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SUMMARY 

A  review  of  recent  work  on  turbine  heat  transfer  performed  at  the  RAE  (Pyestock)  is 
presented.  The  work  covers  the  effects  of  secondary  flows  on  turbine  nozzle  guide  vane 
heat  transfer  with  and  without  film  cooling.  It  is  shown  that  the  heat  load  to  the  platforms 
(cndwalls)  are  significantly  affected  by  the  second^  flow  action.  The  platform  film 
cooling  data  has  ^en  well  coirelated  with  flat  plate  single  row  film  cooling  data  to  within 
±11%.  A  three-dimensional  Navier-Stokes  computational  study  of  the  effects  of  turbine 
inlet  temperature  distortion  on  the  thermo-fluid  mechanics  within  a  rotating  blade  passage  is 
given.  It  is  shown  that  the  temperature  distortion  is  modified  within  the  rotor  blade  and  can 
lead  to  increased  pressure  side  and  over  tip  heat  transfer. 
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MgASgRPjEOTS  ^  OOMPOTATIONS  OF  EXTERNAL  HEAT 
TRANSFER  AND  FlUf  COOLING  IN  TURBINES 


S.P.  Harasgama,  C.D.  Burton  «  K.S.  Ghana 
Royal  Aerospace  Establishment 
Propulsion  Department,  Pyestock, 
Farnborough,  Hampshire,  U.K. 


ABSTRACT 

A  review  of  recent  work  on  turbine  heat 
transfer  performed  at  the  RAE  (Pyestock) 
is  presented.  The  work  covers  the  effects 
of  secondary  flows  on  turbine  nozzle  guide 
vane  heat  transfer  with  and  without  film 
cooling.  It  is  shown  that  the  heat  load  to 
the  platforms  (endwalls)  are  significantly 
affected  by  the  secondary  flow  action.  The 
platform  film  cooling  data  has  been  well 
correlated  with  flat  plate  single  row  film 
cooling  data  to  within  ±  11%  .  A  Three- 
Dimensional  Navier-Stokes  computational 
study  of  the  effects  of  turbine  inlet 
temperature  distortion  on  the  thermo-fluid 
mechanics  within  a  rotating  blade  passage 
is  given.  It  is  shown  that  the  temperature 
distortion  is  modified  within  the  rotor 
blade  and  can  lead  to  increased  pressure 
side  and  over  tip  heat  transfer. 


HOMENCLATURE 

A  Film  cooling  superposition  parameter 

B  Film  cooling  superposition  parameter 

-B/A  Cooling  Effectiveness 

Ct  Vane  tangential  chord 

d  Cooling  hole  diameter 

G  Blowing  rate  (FiDi/rgOg) 

h  Heat  transfer  coefficient 

I  Momentum  flux  ratio  ( r^U j^/FgOg^ ) 

k  Thermal  conductivity  of  air 

M  Mach  number 

Nu  Mussel t  number 

Q  Heat  Flux 

Re  Reynolds  number 

S  Streamwise  dimension 

S2  Pitchwise  dimension 

T  Temperature 

U  Velocity 

X  Axial  dimension 

r  Density 

e  Temp.  Diff.  Ratio  (Tg-Ti ) /(Tg-T„) 

*  Cooling  function  parameter 


cp  Corrected  for  property  variations 

g  Gas  path  value 

i  Injectant  value 

o  Datum  (uncooled)  value 

t  Total  (stagnation)  value 

w  Wall  value 


INTRODUCTION 

Turbine  heat  transfer  research  has  been 
undertaken  at  the  RAE  in  both  experimental 
and  computational  fields.  In  the  main,  the 
work  has  concentrated  on  the  heat  transfer 
resulting  from  the  three-dimensional  (3D) 
nature  of  the  flow  field  within  high 
pressure  turbine  stators  and  rotors.  These 
secondary  flows  have  been  well  catalogued 
in  terms  of  their  aerodynamic  impact''^. 

The  endwall  inlet  boundary  layers  give 
rise  to  two  sets  of  secondary  flows.  The 
passage  vortex  is  caused  by  overturning  of 
the  endwall  t>oundary  layers  within  the 
turbine  passage.  In  addition,  a  horse-shoe 
vortex  is  formed  by  the  stagnation  of  the 
endwall  boundary  layers  at  the  blade 
leading  edge.  Both  mechanisms  cause  a 
migration  of  flow  towards  the  suction  side 
of  the  passage  and  can  lead  to  increased 
losses. 

Heat  transfer  studies  on  the  effects  of 
secondary  flow  have  also  been  reported  by 
several  authors^'  .  It  is  shown  that  the 
horse  shoe  vortex  has  a  major  effect  on 
the  platform  heat  transfer  and  also  on  the 
vane  suction  side.  The  horse  shoe  vortex 
strips  off  the  incoming  boudary  layer  and 
this  gives  rise  to  a  new  boundary  layer 
st®rti.jg  just  downstream  of  the  separation 
line.  This  results  in  enhancing  the  heat 
load  in  the  downstream  region. 

The  heat  transfer  effects  due  to  radial 
variations  in  inlet  temperature  are 
brought  about  by  hot  streaks  generated 
within  the  combustor  and  is  essc’i.ially  an 
unsteady  effect'^.  Early  analytical  work’^ 
has  indicated  that  density  st  .-^^tification 
in  a  rotating  pasage  can  gi*  e  rise  to 
secondary  flows  due  to  Cor..olis  and 
Centripetal  accelarations  affecting  the 
flow.  Recent  two-dimensional  viscous 
numerical  simulations'"’®  show  that  the 
hot  gas  is  transported  within  the  blade 
passage  and  is  deposited  onto  the  blade 
pressure  side. Viscous  3D  computations’'’® 
additionally  show  that  hot  gas  migrates 
radially  leading  to  enhanced  rotor  tip 
heat  transfer . 

This  paper  reviews  the  work  done  at  the 
RAE  on  experimental  and  computational  heat 
transfer  in  turbines.  Detailed  results  of 
heat  transfer  on  a  fully  annular  geosietry 
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of  the  aerofoils  and  platforms  of  turbine 
nozzle  guide  vanes  along  with  flow 
visualisations  which  corroborate  the  heat 
transfer  distributions  are  presented.  Full 
coverage  film  cooling  data  on  the  endwall 
of  a  high  pressure  vane  are  also  presented 
and  correlated  with  axial  pressure 
gradient.  A  computational  study  using  a 
fully  3D  Navier  Stokes  code®  on  the 
effects  of  combustor  generated  temperature 
distortion  within  a  rotating  turbine  blade 
passage  is  also  reviewed. 

TEST  FACILITY 

The  Isentropic  Light  Piston  Cascade  (ILPC) 
as  shown  in  Figure  1  is  derived  from  the 
short  duration  thermc -fluids  facilities 
pioneered  at  Oxford  University’.  The 
facility  operates  by  rapidly  compressing 
air  with  a  light  piston  which  raises  its 
pressure  and  temperature  to  a  required 
level.  Subsequently  a  fast  acting  valve 
is  opened  and  the  air  is  discharged 
over  the  set  of  turbine  vanes/blades.  The 
compression  time  of  the  facility  is  I.Osec 
with  a  run  time  of  0 . 5sec  when  operating 
at  a  gas-to-wall  temperature  ratio  of  1.5. 
Engine  similarity  conditions  of  Reynolds 
number  and  Mach  number  prevail  for  the 
O.Ssec  run  time.  For  film  cooling  studies, 
mixture  of  air  and  carbon-dioxide  (CO2) 
was  used  to  obtain  the  correct 
coolant-to-gas  density  ratio  whilst 
maintaining  the  correct  gas-to-wall 
temperature  ratio. 

04IE 

'  MAIN  IS;.?*  vA;v6S  -  C**' 


The  vanes  under  test  are  manufacj^ured  from 
machineable  glass  ceramic  (MACOR  )  which 
has  a  low  thermal  diffusivity.  Thin  film 
heat  transfer  gauges  are  deposited  onto 
the  surface  of  the  vanes  and  monitor  the 
temperature  during  the  run  time.  The  heat 
transfer  coefficients  are  obtained  by 
solving  the  one-dimensional  transient  heat 
conduction  equation.  A  vane  with  typical 
Instrumentation  and  platform  film  cooling 
holes  is  shown  on  Figure  2.  A  sufficient 


number  of  thin-film  gauges  were  placed  on 
the  vanes  to  give  accurate  measurements  of 
the  heat  transfer  distributions.  In  all 
cases  aerodynamic  data  were  acquired  at 
the  same  locations  as  the  heat  transfer 


measurements.  This  was  achieved  by  having 
static  pressure  tappings  on  the  vanes  and 


Fig  2  Heat  Transfer  S  Cooling  Blade 


total  pressure  probes  upstream  of  the 
turbine  module.  The  inlet  free  stream 
turbulence  intensity  was  set  to  around  6% 
by  a  bar  grid  placed  4 . 5  chords  upstream 
of  the  cascade. 

TORBIMB  STATtMt  HEAT  TRANSFER  RESULTS 

The  heat  transfer  to  the  aerofoils  suction 
surface  is  shown  on  Figure  3.  The  Nusselt 
number  (Nu)  contours  indicate  the  extent 
of  three-dimensionality  of  the  heat  loads. 
For  all  Reynolds  numbers  and  Mach  numbers 
the  secondary  flow  has  transported  low 
energy  endwall  boundary  layer  flow  onto 
the  suction  surface.  These  boundary  layers 
are  relatively  thick  and  thereby  reduce 
the  root  and  tip  heat  transfer  rates.  The 
flow  visualisation  of  Figure  4  tends  to 
corroborate  the  heat  transfer  patterns.  It 
can  be  seen  that  flow  migrating  accross 
the  suction  surface  appears  to  come  from 
the  endwall  boundary  layers. 

Platform  Nusselt  number  contours  are  shown 
on  Figure  5.  The  heat  transfer  patterns 
are  significantly  non-uniform  and  this 
again  is  attributed  to  the  horse-shoe 
vortex  and  the  passage  secondary  flow.  It 
can  be  seen  that  as  the  Reynolds  number  is 
increased  the  high  heat  transfer  region 
migrates  accross  the  endwall  and  moves 
upstream.  A  similar  effect  is  observed  as 
the  Mach  number  is  increased.  In  all  cases 
the  increase  in  Reynolds  number  results  in 
an  increase  of  Nusselt  number  as  expected. 
This  trend  is  the  opposite  for  increasing 
Mach  number .  The  high  heat  transfer  rates 
in  the  downstream  portion  of  the  passage 
are  caused  by  the  horse  shoe  vortex  which 
strips  off  the  incoming  boundary  layer  and 
causes  a  thin  new  ,  high  energy,  boundary 
layer  to  form. 

The  tla^^  visualisation  of  Figure  6  shows 
the  nature  of  the  secondary  flow  on  the 
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TURBIHE  STATOR  FILM  COOLING  RESULTS 

Experimental  data  on  film  cooling  tests 
within  a  fully  annular  turbine  stator  vane 
passage  are  reported.  The  tests  were  done 
on  the  outer  platform  of  the  vane.  These 
tests  were  performed  to  corroborate  a 
design  philosophy  of  placing  cooling  holes 
on  an  iso-Mach  line.  Such  a  design  will 
generate  a  uniform  blowing  rate  accross 
the  whole  passage.  Consider  the  Momentum 
Flux  Ratio  : 

I  =  (ri*Ui2)/(rg*Ug2)  1 

I  =  (Pti-Pi)/(Ptg-Pg)  2 

In  the  turbine  situation  the  free  stream 
and  cooling  total  pressures  are  fixed  by 
compressor  exit  conditions,  as  is  the 
cooling  static  pressure.  Hence  the  only 
variable  is  Pg,  which  is  constant  along  an 
iso-Mach  line.  Placing  the  film  cooling 
holes  along  such  a  line  will  yield  a 
constant  momentum  flux  ratio,  which  for  a 
fixed  density  ratio  yields  a  uniform 
blowing  rate(G). 

For  the  present  tests  the  vane  Reynolds 
number  was  set  to  approximately  2.6x10® 
based  on  exit  conditions  and  tangential 
chord.  The  Exit  Mach  number  was  0.93  and 
the  gas-to-wall  temperature  ratio  <Tg/T»,) 
was  1.30.  The  engine  cooling  design 
required  a  density  ratio  of  1.8  and  a 
temperature  difference  ratio  (8)  of  1.75  . 


In  order  to  achieve  these  conditions  a 
heavy  gas,  carbon-dioxide,  was  mixed 
with  air  to  simulate  density  ratio's  up 
to  2.0  .The  coolant  was  injected  using  a 
separate  reservoir  with  a  fast  acting 
valve  which  fed  a  plenum  just  below  the 
surface  of  the  platforms  being  cooled.  The 
coolant  was  injected  synchronously  with 
'the  onset  of  main  stream  flow. 

The  Nusselt  number  contours  on  the  cooled 
platform  are  shown  on  Figure  8.  These  data 
are  for  increasing  blowing  rate  (G)  from 
0.53  to  1.11  for  a  fixed  6  of  1.62,  also 
shown  is  the  datum  (uncooled)  result.  It 
can  be  seen  that  as  G  is  increased  the 
Nusselt  nuoibers  decrease  significantly. 
Figure  9  shows  a  similar  Nusselt  number 
pattern  for  G  =  1.11  and  6  increasing 
from  0.87  to  1.62.  Here  again  Nu 
decreases  with  increasing  6.  In  both 
figures  it  is  clear  that  the  suction  side 
of  the  platform  is  receiving  more  cooling 
than  the  pressure  side.  This  is  because 
the  secondary  flow  convects  the  coolant 
towards  the  suction  side  which  is 
evident  from  the  flow  visualisation  of 
Figure  6.  Figure  10  shows  the  ratio  of 
cooled- to- uncooled  Nusselt  numbers  for 
G  =  1.11  »  e  =  1.44,  it  is  clear  that  over 
most  of  the  platform  the  heat  load  is 
reduced  by  75%  whilst  at  the  pressure  side 
corner  the  reduction  is  about  20%. 


These  data  have  been  analysed  using  the 
superposition  approach'®  to  film  cooling. 
All  data  have  been  normalised  to  constant 
property  conditions"  using: 


(Nu/Nuo)cp  =  (Nu/Nuq)  *  (T„/Tg)"  3 

The  superposition  representation  is: 
(Nu/Nuo)cp  =  A  +  Be  4 


Where  n  =  -0.25  for  a  turbulent  boundary 
layer.  A  and  B  are  constants  of 
proportionality  in  the  superposition 
formulation' ^ . 


Figure  11  shows  the  linearity  of  the 
energy  equation  when  applied  to  these  data 
and  indicates  that  the  superposition 
formulation  holds  even  in  the  presence  of 
strong  secondary  flows.  This  enables  the 
present  data  to  be  extrapolated  to  other 
engine  operating  conditions.  Also  shown 
are  the  platform  cooling  effectiveness  on 
Figure  12atb.  Figure  12a  shows  the 
streamwise  variation  of  effectiveness 
which  indicates  that  there  is  a  gradual 
reduction  of  -B/A  in  the  downstream 
direction  as  the  coolant  mixes  with  the 
mainstream  gas.  It  also  indicates  a 
reduction  of  -B/A  from  suction  to  pressure 
side  which  is  more  explicitly  presented  in 
Figure  12b.  This  indicates  that  the 
suction  side  effectiveness  is  greater  than 
the  pressure  side  value  for  any  axial 
distance  due  to  secondary  flow  effects. 

The  present  data  have  been  successfully 
correlated  with  single  row  film  cooling 
results  on  a  flat  plate  with  zero  axial 
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pressure  gradient.  This  was  achieved  by 
modifying  the  present  data  with  a 
correction  for  the  turbine  pressure 
gradient.  It  has  been  shown'^  that  the 
ratio  Nu/Nuq  can  be  put  as: 

(Nu/Nuo)cp  =  ♦(S/d)  X  (Uj/Ug)"*/^  5 

A  modification  was  made  to  the  present 
data  by  multiplying  it  by  (U/Uinlet)“' 
and  plotting  against  Equation  5.  The  value 
n=0.2  was  chosen  by  arguing  that,  usually, 
heat  transfer  coefficients  can  be  put  as 
functions  of  Reynolds  number  (velocity) 
raised  to  the  power  1/5th  or  4/5th.  The 
results  of  this  analysis  are  shown  on 
Figure  13aSb.  Figure  13a  shows  the  data 
plotted  against  the  discrete  single  row 
injection  parameter  (Eq-5).  It  is  seen 
that  (Nu/Nuo)cp  lies  below  the  flat  plate 
zero  pressure  gradient  correlation.  The 
corrected  data  are  shown  on  Figure  13b.  It 
can  be  seen  that  Equation  5  now  represents 
the  platform  data  very  well.  The  scatter 
about  the  line  is  only  +  11%,  which  is 
considered  to  be  very  good  in  the  present 
case  of  large  secondary  flows. 


fig  1241  ttreaafte*  vectsttsm  ef  iff  f*<t  !«•«>«••  >  $/c, 
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HOT  GAS  MIGRATION  WITHIN  ROTOR  BLADES 

The  effect  of  inlet  temperature  distortion 
on  the  thermo-fluid  mechanics  in  the  rotor 
blade  has  been  computed''  using  a  steady 
state  three-dimensional  Navier-Sto)tes 
code®.  The  inlet  temperature  distortion  is 
generated  by  the  combustor  burners  and 
affects  the  rotor  blading  as  they  pass 
through  these  hot  patches;  it  is  therefore 
an  unsteady  effect.  However,  an  "average" 
effect  of  the  hot  gas  migration  can  be 
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evaluated  using  a  steady  state  3D  viscous 
computation  in  order  to  gain  an 
understanding  of  the  flow  physics. 

The  simulation  of  inlet  temperature 
distortion  was  achieved  with  a  parabolic 
radial  temperature  profile  as  shown  on 
Figure  14.  The  computational  grid  for  the 
rotor  blade  is  shown  on  Figure  15.  The  3D 
Navier-Stokes  code®  was  run  with  67-axial, 
25-radial  and  25-blade-to-blade  grid 
points.  The  predicted  secondary  flow 


Pig  14  Inlgt  Radial  Tanperature 


patterns  are  given  on  Figure  16.  It  can  be 
seen  that  the  computations  with  radial 
temperature  distortion  give  higher 
secondary  flows.  Figure  17  shows  the  over 
tip  Mach  number  distributions  which  reach 
a  peak  value  of  1.80  and  this  will  also 
enhance  the  rotor  tip  heat  transfer  rates. 


Contours  of  rotor  relative  total 
temperature  are  given  in  Figure  18.  This 
shows  that  the  hot  gas  migrates  from  the 
mid-span  region  at  turbine  inlet  to  the 
presssure  side  tip  of  the  blade  by  80% 
axial  chord.  The  hot  gas  subsequently 
enters  the  tip  gap  and  is  entrained  into 
the  tip  vortex.  The  heat  transfer 
enhancement  on  the  blade  pressure  side  due 
to  hot  gas  migration  is  shown  on  Figure 
19.  This  indicates  a  75%  enhancement  due 
to  hot  gas  migartion.  The  actual  increase 
will  be  periodic  and  fluctuate  between  0% 
and  75%  as  the  blade  passes  through  the 
hot  gas  patches. 

CONCLUSIONS  >  OBSERVATIONS 

Full  three-dimensional  heat  transfer 
distributions  have  been  measured  on 
turbine  vanes  operating  at  engine 
representative  conditions.  The  effects  of 
secondary  flow  on  heat  tranfer  are  large. 
Rotor  blade  heat  transfer  is  significantly 
affected  by  inlet  temperature  distortion. 
The  main  conclusions  are  as  follows: 

-  The  secondary  flow  reduces  the  root  and 
tip  heat  transfer  on  the  vane  suction 
side  due  to  low  energy  endwall  fluid 
migrating  from  the  endwall  regions. 

-  Endwall  heat  transfer  is  also  affected 
by  secondary  flow  leading  to  high  heat 
transfer  near  the  pressure  side  trailing 
edge . 

-  Flow  visualisations  indicate  that  the 
horse-shoe  vortex  sweeps  off  the 
incoming  endwall  boundary  layer. 

This  leads  to  high  energy  free  stream 
gas  forming  a  thin  new  boundary  layer, 
giving  high  heat  transfer. 

-  Platform  film  cooling  data  indicate  good 
effectiveness  when  the  injectant  emerges 
into  a  uniform  pressure  field.  This  data 
has  been  successfully  correlated  to  ±11% 
of  flat  plate  results. 
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Figure  1?  K»ch  KuMbcri  In  Tip  Cap,  Axial  view 


-  Rotor  blade  pressure  side  heat  transfer 
is  increased  due  to  hot  gas  migration. 
The  maximum  increase  can  reach  75%  as 
the  blades  pass  periodically  through  hot 
gas  patches  generated  by  the  combustor. 
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